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Using  high- resolution  time- sequence  photographs  of  solar  granulation  from  the  SOUP  experiment  on 
Spacelab  2,  we  have  observed  large-scale  horizontal  flows  in  the  solar  surface.  Our  measurement  method  is 
based  upon  a  local  spatial  c.-oss  correlation  analysis.  The  horizontal  motions  have  amplitudes  in  the  range 
300  to  1000  m/s.  Radial  outflow  of  granulation  from  a  sunspot  penumbra  into  surrounding  photosphere  is  a 
striking  new  discovery.  We  see  both  the  supergranulation  pattern  and  cellular  structures  having  the  scale  of 
mesogranuJation.  The  vertical  flows  that  are  inferred  by  continuity  of  mass  from  these  observed  horizontal 
flows  have  larger  upflow  amplitudes  in  cell  centers  than  downflow  amplitudes  at  cell  boundaries. 
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LARGE-SCALE  HORIZONTAL  FLOWS  FROM 
SOUP  OBSERVATIONS  OF  SOLAR  GRANULATION 


L.J.  November  (NSO/SP),  G.W.  Simon  (AFGL/SP),  T.D.  Tarbell, 
A.M.  Title,  wid  S.H.  Ferguson  (LPARL) 


ABSTRACT 

Using  high- resolution  time-sequence  photographs  of  solar  granulation  from  the  SOUP  experiment  on 
Spacelab  2,  we  have  observed  large-scale  horisontal  flows  in  the  solar  surface.  Our  measurement  method  is 
based  upon  a  local  spatial  cross  correlation  analysis.  The  horitonUd  motions  have  amplitudes  in  the  range 
300  to  1000  m/s.  Radial  outflow  of  granulation  from  a  sunspot  penumbra  into  surrounding  photo^here  is  a 
striking  new  discovery.  We  see  both  the  supergranulation  pattern  and  cellular  structures  having  the  scale  of 
mesogranulation.  The  vertical  flows  that  are  inferred  by  continuity  of  mass  from  these  observed  horisontal 
flows  have  larger  upflow  amplitudes  in  cell  centers  than  downflow  amplitudes  at  cell  boundaries. 


1.  Introduction 

Since  the  Doppler  effect  measures  line-of-sight 
velocities  it  cannot  be  used  to  determine  horizontal 
motions  at  disk  center.  The  only  way  to  observe 
such  flows  is  through  proper  motion  studies  using 
solar  feattires  as  tracers.  Granules  are  the  obvious 
choice  as  tracers  of  photospheric  flows,  but  because 
of  their  small  size,  coupled  with  the  blurring  and 
geometric  distortion  due  to  atmospheric  tur¬ 
bulence,  they  have  proven  to  be  lifficult  to  use  in 
ground-based  observations  (Simo.-i  1967).  During 
its  10  minute  lifetime,  a  typical  granule  might 
move  only  a  few  hundred  kiioraeiers,  about  one- 
third  of  an  arcsec,  which  is  much  less  than  the  size 
of  the  random  fluctuations  caused  by  atmospheric 
seeing,  even  under  the  best  conditions. 

Now,  however,  with  the  high-quality 
distortion-free  movie  sequences  from  SOUP,  we  can 
meMure  granulation  proper  motions  with  high 
accuracy.  The  SOUP  image  quality  was  so  good 
(pointing  jitter  was  about  0.003  arcsec,  or  2  km  on 
the  sun)  that  we  can  measure  flow  speeds  to  an 
accuracy  of  15  -  20  m/s  nns.  This  is  much  less 
than  errors  in  the  flow  measurements  due  to  ran¬ 
dom  motioiu  and  evolution  of  individual  granules 
and  the  residual  efi'ecta  of  the  large-scale  five- 
minute  oscillations.  Analysis  of  these  various 
"solar  noise"  contributors  indicates  an  overall  error 
of  about  80  m/s  per  independent  spatial  sample  in 
the  results  to  be  discussed  below.  This  is  far 
smaller  than  the  300  -  1000  m/s  signals  which  we 
rqwrt. 

Some  of  the  preliminary  results  have  already 
been  reported  (Title  et.al.  1986,  November  et.al. 
1986b). 

2.  Technique 

Proper  motion  is  determined  from  these  data 
by  computing  the  displacement  that  maximizes  the 
cross-correlation  between  two  image  Gelds 


separated  by  a  time  At.  In  our  analysis  we  have 
used  values  of  At  from  10  to  60  s.  The  images 
are  first  filtered  to  remove  some  photographic  grain 
noise  and  any  large-  scale  intensity  gradients  so 
that  the  remaining  signal  is  mainly  the  solar 
granulation.  Then  the  cross-correlation  function  is 
competed  at  each  "locale"  on  the  images.  The 
locale  around  each  image  point  is  defined  as  the 
product  of  the  intensity  field  and  a  gaussian  win¬ 
dow  centered  on  that  point.  The  gaussian  window 
used  in  this  analysis  had  a  full  width  at  half  max¬ 
imum  (FWHM)  of  4.2  arcsec.  The  vector  displace¬ 
ment  and  transverse  flow  velocity  of  the  features 
within  the  window  during  the  interral  At  are  then 
found  by  locating  the  cross-correlation  maximum. 
The  window  is  moved  from  pixel  to  pixel  over  the 
entire  image  to  give  the  vector  displacement  at 
every  pixel  locale.  This  computational  method  is 
described  by  November  (1986a).  The  4.2  arcsec 
size  of  the  gaussian  was  found  to  be  a  reasonable 
compromise  between  sensitivity  and  resolution  in 
the  final  map  of  transverse  flow  velocitiea  We 
have  developed  a  similar  algorithm  at  Lockheed 
that  uses  the  minimum  of  a  squared  residual  func¬ 
tion.  The  agreement  in  displacements  measured 
by  the  two  techniques  is  excellent. 

Figure  1  is  an  example  of  a  vector  displace¬ 
ment  map  formed  by  averaging  160  original  vector 
maps  spanning  1600s  of  time.  The  window  rise  is 
plotted  in  the  lower  left  of  each  figure;  this 
represents  the  actual  spatial  resolution  of  the  dis¬ 
placement  measurements.  Note  that  we  consider¬ 
ably  oversample  this  resolution  in  the  analysis  and 
display  vectors  formed  from  overlapping  granula¬ 
tion  fields. 

We  have  calibrated  this  technique  by  cross- 
correlating  a  digitized  photographic  image  with 
itself  after  applying  a  Imown  spatial  shift  in  the 
microdensitometer  digitizing  stage.  The  precision 
depends  very  sensitively  upon  the  method  used  for 
interpolation  in  the  cross-correlation.  With  the 
data  sampled  at  or  exceeding  the  diffraction  limit 
of  the  telescope,  a  2-D  quadratic  interpolation  to 
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Figure  1.  Average  proper  motion  vector  map  -  quiet 
region.  The  proper  motion  at  each  spatial  locale  has  been 
defined  by  spatial  cross-correlation  analysis  withm  a  sliding 
gaussian  window.  The  granulation  is  shown  with  super¬ 
posed  vectors  having  lengths  proportional  to  the  amplitude 
of  the  local  Dows.  The  scale  of  the  vectors  is  given  by  the 
arrow  shown  below  these  figures.  The  gaussian  window  size 
of  4.2  arcsecs  FWHM  is  indicated  by  the  size  of  the  circle  in 


Figure  3.  Average  proper  motion  vector  map  -  pore 
region  This  is  a  region  containing  numerous  pores.  Notice 
that  the  amplitude  of  the  flows  is  considerably  reduced  in 
the  vicinity  of  the  pores.  Away  from  the  pores  cellular  flows 
characteristic  of  those  seen  in  the  quiet  sun  image  are  again 
apparent. 
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Figure  4.  A  magnified  portion  of  the  quiet  proper-motion 
map  from  Figure  1  is  magnified  and  more  vectors  are 
displayed.  This  demonstrates  the  cellular  nature  that  is 
characteristic  of  the  large-scale  proper-motion  flows. 


Figure  2.  Average  proper  motion  vector  map  -  sunspot 
region.  A  large  isolated  spot  was  photographed  by  SOUP, 
and  shows  radial  outflow  from  the  penumbra  and  nearby 
photosphere  No  vectors  are  shown  in  the  umbra  of  the 
sunspot  as  there  was  inadequate  signal  there  for  the  cross- 
correlation. 


the  crose-corre^i-.!!  ,  r-Tes  «  meejurement  of  the 
tinege  ihift  i  .  oia;  ezrcr  of  less  th&n  3  per¬ 

cent  for  displa*;e=ic:3  -  —  'Jse  rsa^e  0.005  pixel  (0.6 
km  on  the  sun  ;c  ;  nacei 

The  error  c  ue  Tme*snreiaent  is  estimated  by 
examining  the  erznxzjcn:-.  of  the  rector  displacement 
in  time.  We  corrj'ue  ;  tane-XTeraged  vector  dis¬ 
placement  maf  m  n  rtgnre  1:  the  fluctuations 
about  this  averag*  {rv«  upper  estimate  of  the 
noise  in  the  rsear-iTKEmer.;.  The  Fourier  power 
spectrum  of  vi-»  lanmncral  variation  shows  two 
components:  a  jc»-rwaa«!cy  feature  and  a  smaller 
high-frequency  notse  tall.  The  tail  reflects 

the  statistical  uut.guunrrry  m  the  measured  dis¬ 
placement  and  a  -pnzuxibrr  due  principally  to  pho¬ 
tographic  granu.a.-u-’  -.noise  It  coittributes  an  rms 
error  of  20  m  s  cr  :  1.1.  pixel,  to  the  160  frame 
average  of  v-ectrr  ngpiacements.  The  low. 
frequency  compcnc::  nas  maximum  power  at 
about  200s  period  inu  .  conadbuiea  an  rms  error  of 
80  m/s  to  the  avEnze.  The  level  of  this  "solar 
nobe"  may  vary  quiet  sun  to  plage  or  sun¬ 

spot  regions. 

1.  .-.nnalyiis 


and  from  (34,17)  to  (28,24).  We  can  understand 
this  structure  by  analysing  the  divergence  of  the 
horisontal  flow  n  *  (u,v),  y-n  —  du/dx  +  &v/dy. 
The  divergence  b  shown  in  Figure  7  as  contours 
superposed  on  the  quiet  granulation  field  from  Fig¬ 
ure  1.  Here  we  identify  the  outflow  at  (34,17)  m 
the  large  positive  divergence  that  reaches  an 
amplitude  of  5.510~4  $-1.  Thb  outflow  b  lui^ 
rotmded  by  smaller  cells.  Apparently  the  superpo¬ 
sition  of  the  flows  from  the  smaller  cells  creates  the 
impression  of  the  radial  velocity  structure. 

One  large  sunspot  was  observed  in  the  area 
photographed  by  SOUP.  Figure  2  shows  thb  region 
with  superposed  flow  vectors.  Radial  outflow 
extends  from  the  edge  of  the  sunspot  penumbra 
into  surrounding  photosphere  with  amplitude  300  - 
700  m/s.  Thb  annulus  or  "collar"  around  the 
penumbra  has  a  width  of  about  five  arcsec.  Thb 
outflow  is  reduced  adjacent  to  two  small  sectors  of 
penumbra,  one  of  which  b  poorly-formed  and 
intermediate  in  brightness  between  normal  penum¬ 
bra  and  photosphere  We  show  in  Figure  8  the 
sunspot  region  with  superposed  contours  of  flow 
divergence.  Away  from  the  sunspot  we  find  cellu¬ 
lar  inflows  and  outflows  similar  to  those  seen  in  the 
quiet  region  of  Figure  1. 


Three  spec:.£:  rcrcmcns  of  tie  166x250  arcsec, 
160  frame,  27  inata  set  have  been  analysed: 

a  region  of  qiue;  jrnt.  ar.  isolated  susspot,  and  a 
pore  region.  Thu  tati  .  set  was  obtained  between 
19:10  and  19:2-5  TI  zz  I  Angnst  1985,  in  the 
vicinity  of  NOAA  acctrve  region  4682.  Single 
frames  of  these  imigv  seids  are  shown  with  their 
overlying  averaxe  T»n:— r-  proper  motions,  respec¬ 
tively,  in  Figctres  1  .  and  3.  Each  of  these 

images  was  prccuac  digitbing  piositive  copies 
of  the  original  £.':rr,i  :  us  look  at  these  figures  in 

some  detail. 

In  the  qui«t  sttn  ■  nnciograph  (Figure  1)  we 
identify  two  regiora  :i  larger  horisontal  outflow 
emanating  from  zrLzers  s  a:  (34,17)  and  (8,25)  (that 
b  to  say,  the  coor:  rmaies  of  the  outflow  center 
(34,17)  are  at  x=34  arcasec  and  y=17  arcsec  in  the 
figure).  We  have  nagxnified  the  area  centered  at 
(8,25)  and  show  it  n  r'zirnre  4.  Here  we  see  abo  a 
nearby  linear  scnk  regnnor  around  (15,22).  Flows 
which  emanate  fros  funngle  points  and  then  flow 
toward  and  along  a  lue  ?  are  fairly  characterutic  of 
quiet-  region  veaociry  taiaisems.  However,  note  the 
nice  exception  sriovni  convergence  to  a  point  at 
(25,34)  in  Figure  I  Tm.j  small  region  b  magnified 
in  Figure  5.  The  icaie  of  many  of  the  cellular 
features  b  typicaly  11  arcsec:  we  associate  thb 
with  mesogranu,atiai  ?'--»avemb«r,  et,al.  1981). 

The  region  amrnc.  134.17)  in  Figure  1  b 
shown  enlarged  m  Fifunr”  6.  Thb  feature  has  the 
typical  supergrairuif  sea.*,  of  about  30  -  40  arcsec 
(Simon  and  Leixbiru  .J-&64  The  outflow  in  thb 
example  appears  t:  ju-u-  abrupt  cheuiges  in  flow 
direction  partitirue;  s.ii;ucr  rad.aJ  lines  approxi¬ 
mately  45  degrees  ip-mr  zzit  the  lines  of  flow 
from  (34,17)  tc  2;  irzzz  ,'34.1“)  to  (25,17), 


Figure  3  shows  a  region  containing  numerous 
pores  with  superposed  flow  vectors.  The  granule 
proper  motions  are  much  reduced  in  the  pore 
region,  but  away  from  the  pores  we  again  find  cel- 
lulair  inflows  and  outflows  as  in  the  quiet  sun.  The 
lack  of  flow  in  the  pore  region  b  iltustrated  by  the 
small  amplitude  of  the  divergence  signal  over  the 
pores  shown  in  Figure  9.  The  region  in  the 
immediate  vicinity  of  the  pores  b  magnified  in 
Figure  10,  and  here  the  reduced  magnitude  of  flow 
vectors  can  be  clearly  seen.  The  flow  velocities  in 
the  pore  region  are  about  100  m/s,  about  the  sue 
of  the  "solar  nobe"  in  our  measurements,  and  thus 
diflScult  to  observe. 


4.  Discussion 

In  the  divergence  maps  of  the  cellular  quiet 
sun  Sows,  we  find  that  the  places  of  positive  diver¬ 
gence  characterbtically  have  larger  amplitudes 
than  the  places  of  negative  divergence  even  though 
in  its  definition  the  average  divergence  over  the 
area  of  the  image  b  set  to  sero.  Compare,  for 
example,  the  extremes  of  divergence  in  Figure  1  at 
region  (8,25)  with  5.5- 10-4  s-i  and  the  region  at 
(24,31)  with  divergence  -4.5-10-4  s-i.  We  quan¬ 
tify  the  difference  between  upflow  and  downflow 
regions  by  plotting  a  hbtogram  of  the  divergence. 

It  extends  to  larger  amplitudes  for  the  positive 
divergences.  The  mean  of  the  positive  divergence 
b  1.29-10-4  s-i  and  the  mean  of  the  negative 
divergence  b  -1.17-10-4  s-i;  thus  the  sinks  occupy 
more  area  than  the  souit  es. 

The  sinks  and  sources  are  associated  with 
downflows  and  upflows  through  the  mass  con- 

le^  MMe  Cop) 
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Fipire  5.  Magnification  of  another  portion  of  Figure  1  to 
show  a  sink  point.  More  typically,  however,  sources  are  iso¬ 
lated  points,  while  sinks  lie  along  lines  of  convergence. 
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Figure  6.  Magnification  of  a  third  part  of  Figure  1  show¬ 
ing  a  large  supergranular  cell. 


Figure  7.  Divergence  of  quiet-region  flows.  The  contours 
show  the  amplitude  of  the  local  divergence  of  the  horiiontal 
flows  from  Figure  1.  The  contour  interval  is  110-*  s->. 
Point  source  regions  are  those  with  positive  divergence 
(solid  contours),  while  the  linear  sink  regions  have  negative 
divergence  (dashed  contours).  The  rero  contour  is  dotted. 
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8.  Divergence  of  sunspot-region  flows  shown  in  Fig- 
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Figure  9.  Divtrgence  of  pore-region  flows  from  Figure  3 
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Figure  11.  Propagation  of  long-lived  tracers.  "Corks" 
embedded  Ln  the  flow  field  of  Figure  1  trace  out  the  flow 
pattern  from  their  emergence  in  regions  of  positive  diver¬ 
gence  until  the>’  congregate  in  sinks  and  network  regions 
where  the  divergence  is  negative.  The  nine  frames  proceed 
from  -8  to  -r8  hours  with  the  actual  observation  time  (hour 
0)  shown  as  a  superposed  uniform  grid  of  corks  and  granu¬ 
lation  intensity  in  the  cente  frame. 
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Figure  10.  Pore  region  magnified.  This  shows  the  vectors 
in  the  immediate  vicinity  of  the  pores  from  Figure  3.  The 
amplitude  of  the  flow  is  considerably  reduced  here  relative 
to  the  quiet  sun  flows. 


tinuiCT  equation:  d(fiw)/dx  =  -V  (pn),  where  w  is 
the  Terticai  and  a  the  vector  horizontal  velocity. 
We  rewrite  this  as:  (w/p)  (6p/9z)  +  dw/ijz  =  -y-a. 
The  horizontal  density  fiuctuations  at  any  level 
must  be  small  because  the  mean  density  is  defined 
by  the  local  pressure  which  is  mainly  hydrostatic, 
and  by  the  temperature  which  is  approximately 
fixed  at  a  fiven  depth.  Replacing  1/p.dp/dz  with 
1  h^  and  9w/ds  with  w/h»,  where  h^  and  h,  are 
the  density  and  velocity  scale  heights,  respectively, 
we  obtain  w(  l/h,  l/h.  )  =  -y  a.  Letting  (l/h„ 
1  h«)  =  1/hm  where  ha  is  the  mass  flux  scale 
height,  then  assuming  h,  <<  h,  (Schmidt,  et.al. 
1965).  and  using  a  value  of  hp=200  km  for  the 
subphotosphere,  ha  ^  200  km  and  we  obtain 
w=26m  s  as  the  average  of  the  sources  and 
w— -23m  s  as  the  average  of  the  sinks.  These 
values,  smaller  than  the  noise  in  etuiier  Doppler 
measurements,  explain  why  it  has  been  so  difficult 
to  detect  vertictd  velocities  in  supergianulation. 

The  sesise  of  our  flow  fields  is  like  that  of  typi¬ 
cal  Benard  cells  that  have  greater  central  upflows 
than  downflows  at  the  boundaries  (Chandrasekhar 
1961).  Further,  the  small  amplitudes  which  we 
infer  are  in  general  agreement  with  Doppler  meas¬ 
urements  of  the  vertical  flows  in  the  supergranula- 
tion.  Simon  and  Leighton  (1964)  were  unable  to 
detect  vertical  motions  at  disk  center  although 
theu'  sesisitivity  was  better  than  100  m/s.  More 
s«iaitive  Doppler  measurements  (Frazier  1970, 
Worden  and  Simon  1976,  Giovanelli  1980, 
November,  et.al.  1982,  Kuveler  1983)  give  vertical 
amplitudes  similar  to  those  we  infer  here.  Our 
proper  motion  measurements  are,  however,  in 
qualitadve  disagreement  with  the  Doppler  observa¬ 
tions  which  showed  larger  amplitude  downflows 
than  upflows.  Thus  they  probably  reflect  another 
effect.  The  difference  might  be  explained  as  a  sys¬ 
tematic  error  in  the  Doppler  signal  caused  by  a 
possible  contrast  effect  in  granulation  over  mag¬ 
netic  elements  (Giovanelli  1980).  More  likely,  the 
Doppler  observation  could  result  from  a  superposi¬ 
tion  of  photospheiic  flows  that  are  dominated  by 
magTietic  rather  than  convective  processes  in  the 
network  boundaries  (Giovanelli  1980,  Gebbie, 
eual.  1981). 

The  large-scale  flow  fields,  both  as  observed 
visually  from  the  movies,  and  quantitatively  from 
the  calculated  vector  displacements,  have  the 
proper  size  and  appearance  of  super-  and  meso- 
granular  convection  Bows.  Therefore  we  would 
expect  these  flows,  at  subsurface  levels  in  the  sun, 
to  move  magnetic  fields  to  the  cell  boundaries,  and 
form  the  loci  of  the  photospheiic  and  chromos¬ 
pheric  networks  The  SOUP  white-light  movies 
indeed  give  us  the  definite  impression  that  Iwes  of 
brigh  ter- than  average  granules  occur  at  these 
places  (converging  sinks  and  line  flows)  where  the 
divergeice  is  negative,  i.e.  at  the  cell  boundaries 
This  bright  photospheiic  network  always  seems  to 
surrour.d  the  upwellings,  or  boiling  centers,  which 
mark  the  darker-appeaiing  super-  and  mesogranule 
cmters  Unfortunately  we  have  not  yet  been  able 
to  locate  from  ground-based  data  any  good  quzdity 


magnetogram  or  chromospheric  spectroheliogram 
coincident  with  these  SOUP  data.  These  are 
needed  to  make  a  definitive  comparison  with  our 
flow  fields.  SOUP  data  from  earlier  and  later 
orbits  are  currently  being  reduced.  Some  of  these 
do  coincide  with  National  Solar  Observatory 
images  of  sufficiently  high  quality  to  make  useful 
comparisons  We  hope  to  analyze  these  data  in 
the  next  few  months. 

Implicit  in  our  analysis  is  the  assumption  that 
the  granules  float,  like  corics,  on  the  supcrgranular 
flows,  and  thus  can  serve  as  tracers  to  define  the 
larger  structures.  Suppose,  then,  that  we  start 
with  a  uniform  distribution  of  granule  "corks" 
superposed  on  and  subject  to  the  pattern  of  arrows 
which  define  our  average  flow  fields.  As  time 
evolves,  the  corks  move,  and  as  they  reach  the  next 
pixel  point,  they  are  subjected  to  the  slightly 
different  flow  field  corresponding  to  that  location. 
We  have  made  a  movie  of  such  cork  flows  in  which 
we  extrapolate  our  measured  velocities  forward  and 
backward  in  time  assuming  that  they  do  not 
change  over  the  36  hour  length  of  the  movie.  It  is 
interesting  to  note  that  the  corks  form  into  a  net¬ 
work  pattern  in  two  to  three  hours;  some  of  the 
corks,  as  expected,  congregate  at  the  few  sink  holes 
in  our  images.  If  we  run  such  a  cork  movie  in 
reverse,  the  corks  quickly  concentrate  in  the  upwel¬ 
lings  which  form  the  meso-  and  supergranule 
centera  We  illustrate  the  cork  motions  in  Figure 
11  where  we  show  the  movie  frames  from  times 
-8,  -6,  -4,  -2,  0,  2,  4,  6,  and  8  hours  relative  to 
the  actual  time  of  observation  for  the  quiet  sun 
image  of  Figure  1.  WTiile  these  cork  movies  give 
only  an  artificial  picture  of  the  development  of 
large-scale  flows,  they  do  agree  well  with  our  intui¬ 
tion.  Fortunately  there  exist,  in  all,  about  16 
hotirs  of  SOUP  data  so  that  it  will  be  possible  to 
make  an  authoitic  study  of  the  evolution  of  these 
flows. 
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